Hamiltonian theory of the half-filled Landau level with disorder: 

Application to recent NMR data 

Ganpathy Murthjc 
Department of Physics and Astronomy, 
University of Kentucky, Lexington KY 40506-0055 

R. ShankaiH 

Department of Physics, Yale University, New Haven, CT 06520 

Abstract 

The Hamiltonian Theory of the fractional quantum Hall effect is an operator description that 
subsumes many properties of Composite Fermions, applies to gapped and gapless cases, and has 
been found to provide results in quantitative accord with data on gaps, relaxation rates and po- 
larizations at temperatures of 300mK and above. The only free parameter is A, which is related 
to the sample thickness and appears in the Zhang-Das Sarma potential v(q) = ^^-e~ qlx where I 
and k are the magnetic length and dielectric constant. Here we examine the recent data of Tracy 
and Eisenstein on the nuclear magnetic resonance relaxation rate at filling factor v = \ deduced 
from resistivity measurements at temperatures as low as 45mK. We find that their results can be 
satisfactorily described by this theory, if in addition to a v(q) with A ~ 2, a constant disorder width 
T ~ WOmK is incorporated. 
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I. INTRODUCTION 

There are several theoretical approaches to the fractional quantum Hall effect (FQHE) [l[ 
For fractions of the form v = one can write down inspired trial wave functions 
following Laughlin^ to describe the incompressible liquid ground state and its gapped 
excitations. These excitations can be further studied using an extension of the Bijl-Feynman 



approach of Girvin, MacDonald and Platzman 



4]- 



For a more general class of fractions with filling v = pj (2ps + 1), one can follow Jain 5, f| 
and work with Composite Fermions (CF) which are obtained from electrons by attaching 2s 
flux quanta 3, S|- On average, these attached flux quanta neutralize enough of the external 
flux so that the total is just right for the CF's to fill exactly p Landau levels. When the 
wave functions of these unique states are gauge transformed back to electronic language 
and projected to the Lowest Landau Level (LLL) of electrons, they become excellent trial 
wave functions. The quasiparticle excitations can likewise be described by starting with 
particle-hole excitations of these p-filled CF-LL's (Composite Fermion Landau levels). 

For the Laughlin fractions Zhang, Hansson, and Kivelson Q| provided a microscopic 
Chern-Simons (CS) field theory directly linked to the original Hamiltonian. They converted 
electrons to bosons by adding 2p+l flux quanta so that on average the bosons saw no flux. 
Lopez and Fradkin 10|| extended the CS field theory to the Jain fractions p/(2ps + 1). In 
both the bosonic and fermionic cases fluctuations on top of the mean field description were 
described a CS gauge field. The path integral approach allowed one to more readily compute 
time and space dependent response functions and go to nonzero temperature T. 

The CS theory applied to v = \ leads to CF's that see, on average, zero effective field. 



Kalmeyer and Zhang 



and Halperin, Lee, and Read (HLR)[12| studied the Fermi sea 



of CF's that coupled to CS gauge field. The latter authors performed an exhaustive and 
definitive study making predictions in agreement with many experiments. 

However all the CS theories have one weakness: They do not yield a smooth limit as 
one projects to the LLL by sending m — > 0. This sends the electronic cyclotron frequency 
uj c = — — ► 00, thereby forcing all the electrons to be in the LLL (for v < 1). We know 
that a smooth limit must exist, and that in this limit the kinetic energy should become 
an ignorable constant and the entire Hamiltonian should be just the Coulomb interaction 
v(q) between electrons. The quasiparticle mass, its band structure, its residual interactions, 
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response functions, and so on should be determined solely by v(q). In the Chern-Simons 
theories sending m to zero causes unavoidable singularities. These theories also lead at the 
mean-field level to a CF of charge e and not the correct quasiparticle charge of 

e n 



2ps + 1 

where 



e(l - <?) (1) 



2ps 

c = — - — 2 

2ps + 1 v ; 

is a constant that will appear repeatedly in this paper. 

The electric charge does not come out right because the CS procedure attaches the flux 
but not the screening charge associated with the nucleation of a 2s-fold vortex so evident in 
the wave function description [3], Q]. (In other words, the flux captures only the phase of the 
vortices but not their modulus, with its 2s-fold zero.) In particular, at v = 1/2, c 2 = 1, it 



id- 



does not yield a neutral fermion or describe its dipole structure 

J n 

The present authors addressed these problems by developing a Hamiltonian approach|l4 
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17fl in which the LLL limit can be taken naturally, and many properties of the CF, such 
as its reduced charge e* = l/2ps + 1, and its dipolar nature are encoded unambiguously in 
the operator structure. Relegating the details of this theory to the next section, we merely 
point out that it not only reproduces results found by Monte Carlo calculation based on 

n 

trial wave functions [18|, (for potentials that are not too singular at short distances), it also 
furnishes results in quantitative accord with data on gaps, relaxation rates and polarizations 
at temperatures of 300mK and above for gapped and gapless fractions S|. 

The sense in which the Hamiltonian theory accounts for the data needs some 
elaboration Consider the v — 1/2 data of Dementyev et fl^2ll|. and Khandelwal 22], 
who measured the polarization and the NMR relaxation rate for relatively high temperatures 
(down to 300mK). One can easily find a simple CF theory with an effective mass, and if 
needed, a Stoner interaction J, to fit the data in a limited range. There are two shortcomings 
in this approach. First, the origin of the CF mass and interaction are viewed as unrelated 
when in fact they both originate from the electrostatic interaction between electrons. Next, 



as emphasized by Dementyev et a/|2l| and Khandelwal 22|, if one considers all of their data 
globally, for both polarization and relaxation, the values of (m, J) needed for one observable 
are totally incompatible with those mandated by the other. This is because no single pair 
(m, J) can describe the underlying physics. 



On the other hand, in the Hamiltonian theory, the only free parameter is A, a parameter 
related to the sample thickness in the Zhang-Das Sarma potential ^23 1 

v{q) = (3) 

where I = \J~^ is the magnetic length, and k is the dielectric constant. The parameter 
A is extracted from one data point, after which the theory does quite well at predicting 
the rest of the dataflj]]. As we review the theory, it will become apparent how the CF 
band structure and interactions (whose nonstandard functional forms are fully determined 
by v (q)) w iH arise in a correlated way, and vary with the chemical potential, temperature 
etc., to produce varying effective mass and interaction parameters. 

The current work was triggered by recent experiments of Tracy and Eisenstein (TE)!^. 
They measured the NMR relaxation rate ^- at much lower temperatures than previously 
(going as low as 45 mK) and for a wide range of B fields while keeping the filling factor 
at v = |. The relaxation rate was measured by disturbing the system magnetically, and 
measuring the relaxation of the resistance to its equilibrium value. The precise causal 
relationship between the non-equilibrium value of polarization and the non-equilibrium value 
of the resistance is obviously not needed to extract their common relaxation rate. 

We do not have much to say about the region of very high B where the polarization 
has saturated and the spin up and down bands have separated. Here we simply expect the 
following asymptotic behavior: 

T,{B,T)^e^ Ez ~ E * z (4) 

where Ez is the Zeeman coupling, and E* z is the critical Zeeman coupling at which the 
system becomes fully polarized. 

Our focus is on the region before saturation where the two spin bands overlap. What 
does a naive model of noninteracting CF's predict here? One can easily argue that: 

^^K0)| 4 m* 2 ^5 5 / 3 (5) 

where m* is the effective mass (and controls the density of states) and where u(0) measures 
the height of the wave-function in the transverse direction. Two factors contribute to the 
final power of B: (i) The density of states per unit volume for each spin species scales as 
m* ~ B 1 ! 2 (which in turn follows from setting k 2 F /m* ~ e 2 //) (ii) The factor |w(0)| 4 , scales as 



-B 2 / 3 because the thickness t of the 2DEG scales as n ~ B 1 / 3 [25| and |«(0)| 2 integrates 



to unity in the transverse coordinate. Thus one expects that at fixed T, 

Ti ~ B~ 5/3 . (6) 

However the TE data show a nearly flat dependence on B till saturation begins to set in, in 
contradiction with naive CF theory. 

In the Hamiltonian theory the B dependence of T\ is much more complicated. First the 
density of states (and hence m*) varies with the chemical potential since the bands are not 
quadratic. Next the self-consistent band structure depends on another energy scale, the 
Zeeman energy, and on the thickness parameter A, which scales with B as t/l ~ B 1 ^. All 
of this precludes any naive scaling. One simply has to turn the crank to see how things 
depend on B. We did that and here are the main results. 

We found that while a suitable A ~ 2 (at 4 Tesla) could produce acceptable graphs, the 
data clearly called for a nonzero disorder-induced width Y of order lOOmK for all momentum 
states. The existence of such a T and its order of magnitude both are implicit in TE data 
for the equilibrium resistivity^], reproduced in Fig. (pQ). Assuming a Drude form for the 
conductivity 

° = —T 7 

m* 

with e being the electronic charge, and m* being the effective mass of the CF's, one finds that 
the disorder width T = ^ is a fraction of a Kelvin. Assuming that T is roughly independent 
of n and therefore B, fits the data quite well, as seen in Fig. (DQ). In past comparisons to 
data at 300mK and higher jjjj], disorder played a relatively minor role presumably because 
the temperature was much larger than the T of those samples. 

In this paper our goal is limited to understanding the effect of disorder on the NMR 
relaxation rate, and we merely take Eq. ([7]) as a crude estimate of disorder broadening and 



do not attempt to provide a complete theory of 



ongitudinal conductance at v = 1/2, a 
subject that is yet to be fully developed. (See [12j for a treatment of the charged Chern- 
Simons fermion coupled to a gauge field.) 

Figure [2] shows a comparison of our theory to the T\ data at 45mK and lOOmK. We find 
that the numerical values A = 2.00 and V = 100m.fr produce the best overall agreement. A 
5-independent parallel relaxation channel, with rate 5 x 10 _4 s _1 , (representing nuclear spin 
diffusion out of the 2DEG, and estimated from the TE data), has been added to our theory, 
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Figure 1: The data from Tracy and Eisenstein, Ref. j^J], fitted to a l/B curve. The good fit implies 
that the transport time is independent of B or n. 

and our prediction has been normalized to coincide with the data at B = 4 Tesla. Finally 
(B) was scaled from its value at 4 Tesla. 

Disorder seems to mimic the effects of nonzero temperature in many ways. We see a clear 
flattening of the results for T\ below E* z . Above E z , since the bottom of the minority spin 
band is disorder-broadened, there is always some density of states of both spin species near 
the chemical potential, and therefore T\ rises much slower than the activated form of Eq. 

®- 

Figure [3] shows the comparison of theory versus experiment for 1/T\ versus T for several 
B values. As at low T, the theory shows more 5-dependence than is measured, and does 
not do too well past saturation. 

The plan of this paper is as follows. In Section II we give a brief review of the Hamiltonian 
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Figure 2: Comparison of theory and experiment for T\ vs B at A5mK (solid line, circles) ) and 
WOmK (dotted line, squares). The optimal values A = 2 (at 4T) and T = WOmK were employed. 
In the comparison, an "offset" of 5 x lO^s" 1 has been added to the nuclear relaxation rate to 
account for the diffusion of nuclear spins out 2DEG into the bulk. 

Theory. In Section III we introduce our constant width disorder model, set up the HF 
equations, and present a formula for the NMR relaxation rate and polarization. Section IV 
is devoted to a comparison of theory to experiment and the role of our two free parameters, 
and of the temperature, on the NMR relaxation rate. In Section V we address the effects of 
Landau-level mixing qualitatively, and in Section VI we conclude by presenting some of the 
limitations of our work, how they may be overcome, and some open questions. 



II. REVIEW OF THE EXTENDED HAMILTONIAN THEORY 

We shall furnish only a brief summary of the Hamiltonian theory since complete de- 
tails may be found in our review article[17j. Our starting point will be the full electronic 
Hamiltonian 

« = ES + /(05'«)M-q) (8) 

where m is the band mass of the electrons, II e = p e + eA is the velocity operator of 
the electrons of charge — e, v(q) is the Coulomb interaction between the electrons, and 



p e (q) = J2i ex P ~ ' r ei is the electronic density operator. 

To facilitate the projection to the LLL, we diagonalize the kinetic energy by decomposing 
the electron coordinates and momenta into cyclotron (77J and guiding center (R e ) variables 

r) e = i 2 z x n e (9) 

R e = r e - l 2 z x n e (10) 

where / = \Jhc/eB is the magnetic length. These coordinates obey the commutation 
relations 

[flex , ley] = ^ l^ex , Re,] = -U* [r] e , R e ] = (11) 

The kinetic energy depends only on the cyclotron coordinate, with the LLL corresponding 
to the harmonic oscillator ground state for this variable. Expressing the electron coordinate 
as r e = R e + r] e we can now take the LLL limit of the Hamiltonian to get 

H e = J ^v{q)pM)Pe{-^) (12) 
where p e is the electron density operator projected to the LLL. 

p e ( q ) = e -lW 2 ^ e -^. ( 13 ) 
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Figure 3: Comparison of theory and experiment for \/T\ vs T In the comparison, an "offset" of 
5 x 10~ 4 s _1 has been added to the nuclear relaxation rate to account for the diffusion of nuclear 
spins out 2DEG into the bulk semiconductor. The parameters used were A = 2, V = 100m.fr. 
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While we can easily project the Hamiltonian to the LLL, working with it is another thing. 
Here, briefly, are the problems we face and their resolution. The original electronic problem 
assigns to each particle two coordinates (x e , y e ) = r e and two conjugate momenta (p e x,Pey) = 
p e . In the LLL, kinetic energy is quenched but x e and y e become conjugate variables so that 
H = v(x e , y e ) poses a nontrivial quantum problem. 

Equivalently, each electron has a cyclotron coordinate rj e and guiding center coordinate 
R e . In the LLL the Hamiltonian is v (R e ) with R ex and R ey being conjugate variables. 



Quantizing such a problem is tricky and involves working with analytic wave functions 26]. 

Our trick was to begin with new particles, the CF's, equal in number to the electrons, 
each with coordinate r and velocity II = p + e*A, where e* = e/(2ps + 1) is the CF charge. 
Thus II describes particles which will fill exactly p LL's. In what follows we will also work 
with 77 and R, the corresponding cyclotron and guiding center coordinates of the CF. Note 
that the CF variables do not carry subscripts. 

In this space we form the following entities: 

R e = R + rjc = r- ^z x n (14) 
R,, = R + 77/c = r + x n (15) 

The first vector R e is identified with the electronic guiding center coordinate since it obeys 

[R ex ,R ey ] = -il 2 . (16) 

Thus the LLL projected electronic Hamiltonian we set out to solve is just H = v (R e ). Since 
a quantum problem is characterized by the commutation rules, this is a faithful transcription 
of our original mission. However, there is no assurance the degeneracy of our levels will be 
the same now. 

Indeed we can see there is going to be a huge degeneracy because of the other coordinate 
R„ It describes a particle we call the pseudo- vortex. It commutes with R e : 

[Re , = 0. (17) 

and obeys 

[R vx , R vy ] = il 2 /c\ (18) 

Thus H v is a cyclic variable whose dynamics is unrelated to the original problem. Let us 
get acquainted with it anyway, since all this will change. 
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First of all, it follows from its commutation rules that R„ describes a particle of charge 
— c 2 in electronic units. If this object paired with an electron, it would yield an entity with 
total charge e(l — c 2 ) = e/(2ps + 1) = e*. This is exactly how the vortices in the trial wave 
functions screen the electron to give rise to the CfQ|. However the adjective "pseudo" is 
appended for two reasons: (i) The vortices in the wave functions are not creatures with 
their own coordinates independent of electrons, (ii) The pseudo-vortex so far has nothing 
to do with the original problem, but is rather something we introduced when we enlarged 
the Hilbert from the LLL projected case described by just R e to a regular fermionic space, 
i.e., the CF space. 

Since H = v (R e ) does not depend on H v , the dynamics of H v is unspecified. They are like 
gauge variables and their dynamics is determined by gauge fixing. We made the reasonable 
choice that the density p v formed out of H v annihilated all physical states. 

This gives us two options. 

The first is to start with a self-consistent Hartree-Fock (HF) solution to the Hamiltonian 
of Eq. (1T21) written in CF coordinates, and to compute response functions in a "conserving" 



approximation 



28j, such as time-dependent Hartree-Fock (TDHF). In this 



the constraints are satisfied are satisfied "weakly", in correlation functions 



approximation 



29, 



30L 



use 



The second approximation, which we have employed in the pastflHuJQ- 

and 

in this manuscript, is more unconventional. We argue that when acting on exact physical 
states (which are annihilated by p v ) there should be no difference between the projected 
electron density operator p e , and the following preferred density operator 

P P (q) = Pe(q) -c 2 p v (q). (19) 

The new operator p p has the advantage of exhibiting many of the nonperturbative properties 
of the CF without any computation, at the tree level. First, it describes a particle of charge 
e* = e(l — c 2 ) as can be seen by looking at the zeroth order term in an expansion p p (q) in 



powers of q. The order q term 
dictated by Kohn's theorem 



las the correct coefficient to impose another crucial property 
271 ]. that any intra-LLL matrix element of the density operator 
should vanish faster than linear in q as q — > 0. This property emerges from the conserving 
calculation only after some effort It is very striking that one and the same admixture 

of p v , with coefficient c 2 , serves two purposes at once. 

As mentioned this situation is unusual. Normally we have a fixed Hamiltonian and go 
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in search of a trial HF state from a family, here we have a fixed HF state (with p filled LL 
of CF's) and go searching among a set of Hamiltonians equivalent in the constrained space. 
If the constraint is solved exactly, there is nothing to choose between the original p e , the 
preferred p p , or indeed, any arbitrary combination of p e (q) and p„(q). On the other hand in 
HF, the constraint is not satisfied and the preferred combination emerges as the best to use 
since it encodes certain important nonperturbative properties of the CF at leading order. 

We will henceforth use the preferred form of the density in all our calculations, to which 
we now turn. 



III. NMR RELAXATION RATE FOR v = \ 

The assumption that all the nuclei are in thermal equilibrium with each other on a fast 
time-scale leads to the standard Korringa Law for nuclei in contact with degenerate Fermi 
gases|31]. The formalism for the clean v = | system has been previously worked out by one 
of us 19]. We will proceed directly to the case with disorder. Our model Hamiltonian is the 
following: 

H = J -^v(q)p p (ci)p p (-q) - E Z S Z + H dis (20) 

where Ez is the Zeeman energy, H^ s is the coupling to disorder, and we explicitly present 
our preferred density operator for v = \ 

p p (q) = e"^ 2 / 4 "2i sin(^^!)dt s(k ) ds(k ) (21) 

k 

Here d^ s (k) creates a CF in a state with momentum k and spin-projection s. For our simple 
model, we do not need to specify the detailed form of the coupling to disorder. We assume 
that the net effect of H^s in a disorder-averaged treatment is to provide a momentum- and 
spin-independent width Y to every single-particle state. We incorporate this into our HF 
energy calculation as follows. Let the disorder-averaged energy of a state labelled by k, s be 
e s (k) 3 and its occupation be n s (k). Then decomposing the Hamiltonian of Eq. (1201) in the 
HF aproximation, we find 

e,(k) = -E z % + J |fK?)e- (?,)2/2 [l - cos(q x k/ 2 )][l - n s (k + q) - n s (k - q)] (22) 

oo 

n s(k) = / ^e' w ( e /2 +r 2) 1+e ^ (es ( k )_ £ '_^) (23) 
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where (3 = is the inverse temperature. The form of the first equation is the same 

as in the clean system, while the disorder-broadening T makes itself felt in the occupations. 
It says that since a particle in a state of momentum k (defined after averaging) can be in 
a state of energy different from e(k) by an amount of order T, a band of energy states of 
width T will contribute to n(k) by convolution. 

The spectral function can be related to the disorder average of the exact disorder eigen- 
states in the HF basis. Let the exact disorder eigenstates be labelled by a. We can expand 
the operators which annihiliate a CF in a momentum state in terms of a: 

4(k) = ^0 a (k)4« (24) 

a 

Now the retarded single-CF Green's function can be defined as 

Gf(k,k',t) = -iQ(t)(Q\{d s (k,t),d^(k',0)\Q) (25) 
Fourier transforming and expanding in the exact disorder states we find 

Gf(k , k>HE A« (26) 

a 

Now taking the disorder average and comparing to the spectral function, we infer that 

T/tt 



fe(k)fe(ko = 4. r2 + (E ;_ Mk))2 < 27 » 

When Eqs. (l23l) are iterated to self-consistency they provide the HF state of the CF-Fermi 
sea in the presence of disorder. Of course, one needs to maintain the half-filling condition 
at every HF iteration: 

(g^rW + ^MH^ (28) 

Now one needs to address the spin correlations, which are important for the NMR re 
laxation rate 



19(| . The hyperfine interaction between the electrons and the nuclei is given 



by 



y 7e7n ^2|n(0)| 2 I-S e (0) (29) 



where 7 e n are the magnetic moments of the electron and the nucleus, u(0) is the electronic 
wavefunction at the nucleus, I is the nuclear spin, and S e is the electronic spin density 
operator at the position of the nucleus. We define this operator as follows: 

S«(r = 0) = ^J2 e ' ifiXkl2/2 ~ iqir/Ad ^( k - q)^^'( k ) (3°) 



L 2 

k,q 
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Note that this operator does not commute with the projected density operator, but has LLL 
spin-charge commutation relations[32|. Now a Fermi golden rule calculation produces the 
following expression for the NMR relaxation ratejlj]]: 



1 / 87T7e7 n 2 



K0)| 4 | cft<n|5:(r = ) 0)5:(r = ) *)+5*(r = ) 0)52f(r = J 0|n> (31) 



In our evaluation of the spin-correlation function, we will ignore vertex corrections, both due 
to the gauge field (arising from the imposition of the constraints {3, H]) and from disorder 
averaging. We have found from our previous work that the vertex corrections due to the 
gauge field manifest themselves at extremely low energies, where they give rise to a collective 
overdamped mode, but are unimportant at not too low temperatures. The vertex corrections 
due to disorder are important primarily at small q and at energies much smaller than the 
disorder width T. However, our correlator involves a sum over all q, and the experiments 
are carried out at a minimum temperature of AbmK, which is of the same order as T. Thus, 
we expect the effect of disorder vertex corrections to be negligible as well. 

The rest of the calculation is straightforward. Calling the constant 7r( 87ry 3 e7n ) 2 = D and 
performing the disorder average, we have 

J_ = p K°)I 4 v e 4( k ^K+^K)- l i((^-K) 2 +(^-K) 2 ) 

rp Li - 



X 



/ dt[<n|dt T (k 1 ,t)d T (k^,o)|n>(n|d i (k' 1) t)dt i (k 2 ,o)|n> 

— oo 

+ (fi|^(k 1 ,t)d i (k^,o)|n>(n|d T (k' 1 ,t)dt T (k 2 ,o)|n>] (32) 

Expressing the single-particle Green's functions in terms of the spectral function, and inte- 
grating over t, we obtain 

-= tfHgL E I dEn F (E)(l-n F (E))U) x (33) 

J l kika-oo v ' 

+■ T-i ) (34) 



1 

Now we convert the momentum sums into energy integrals by introducing the densities of 
states per unit volume Pn(e) of the | and | spins to get the final expression 

1 = D\u(0)\' J dEdE^e-^^O 10(^1^2 x 
J- 1 



(T*+(E-E r )*)(r*+(E-El)*) [60 > 
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Figure 4: A comparison of the T\ vs B curves for A = 2.00 for different amounts of disorder. The 
dashed line represents the naive CF prediction of B~ 5 ^ 3 below full polarization. T\ is in arbitrary 
units and the curves have been normalized to be equal at 4 Tesla for ease of comparison. 

Here k s = k s (E s ) are to be understood as the result of inverting the energy-momentum 
dispersion relation, and Jo is the modified Bessel function arising from the angular average 
of exp(k| ■ k|/ 2 /2). It is important to note that while the n F (E)(l — n F (E)) factor gives the 
dominant T dependence at low T, there are hidden T dependencies in the energy-momentum 
dispersion relations obtained from the self-consistent HF solution. 



IV. RESULTS AND DISCUSSION 

Let us begin by asking how we extract the values for and T. Note that 1/Ti has the 
form 

± = \u(0)\ 4 f(T,(B),T) (36) 

Let us look at the T dependence of jr at fixed B, say at 4 Tesla. A rough estimate of 
comes from 1 where the following formula is is derived 

P = .13y/B(T) 7/ * u = - P<1 (37) 
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relating the polarization to B(T) in Tesla prior to saturation. From TE data which tells us 
P = 1 occurs first for B{T) ~ 6.25 we deduce = 1.9 as a good start which must be further 
improved by looking at the data. The width Y = lOOmK is found by considering T% versus 
B data for T = 45 mK, T = WOmK. None of these numbers is unique in that the fit is 
not perfect and changing the numbers could improve things in one region and worsen it in 
another. With these values of T and we pick the overall scale, including |w(0)| 4 , to agree 
with the data at 4-Tesla, 45 mK, with the offset of 5 ■ 10~ 4 s -1 from spin diffusion into the 
bulk included. We are now set to make predictions at any other B since T is assumed to be 
constant and the scaling of and |w(0)| 4 are known: 



Fig. [Ushows the effect of T when A is held fixed at 2.00 and the temperature is held constant 
at 45 mK. The dashed line shows the expectations from naive CF theory, which predicts a 
i?~ 5 / 3 dependence until full polarization, and an exponential increase in 7\ thereafter. As 
can be seen, even without disorder, the predictions of self-consistent HF vary substantially 
from the naive CF expectation. However, the qualitative features of the naive expectation 
are indeed present in the clean limit. 

The effect of disorder is to reduce the dependence on B below full polarization, and to 
make the rise after full polarization gentler than exponential. The second feature can be 
understood by observing that even when the disorder-averaged position of the bottom of 
the minority spin band is above the chemical potential, the Lorentzian broadening of the 
band bottom will produce some density of minority spin states at the chemical potential. 
Assuming a constant Y this DOS can be seen to go as (B — £> c ) _2 , which implies a behavior 



In Fig. [5] we show the behavior of 7\ in the clean system (r = 0) at three different 
temperatures. It is clear that at higher temperatures, T\ is lower due to a greater part of the 
band being partially occupied. However, the rise of Ti beyond full polarization is extremely 
rapid, and it is clear that this is not compatible with the experimental dataf^. 




(38) 



(B - B c f for T] 
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A. Polarization 



As mentioned earlier, a compelling reason for preferring the Hamiltonian theory over 
more phenomenological CF theories without a direct link to the electronic Hamiltonian arises 
when one simultaneously tries to explain data on independent quantities like polarization 
and 1/Ti. Where the simpler models call for mutually incompatible values of mass and 
interaction parameters, the Hamiltonian theory, given just A from one data point, is able to 
account for all the salient features of the data. 

To this end we are presenting here the predictions of our theory on polarization P(B,T) 
for the optimal values A = 2, T = lOOmK. Since the dependence of P on T is very weak for 
the range considered by TE, we present just the B dependence in Fig. [61 Unfortunately TE 
do not measure the polarization. We present these results in the hope that one day they 
may be tested. 




- T=50mK 
-T=75mK - 
T=100mK - 



B (Tesla) 



Figure 5: The predictions of the HF approximation of the extended Hamiltonian theory for the 
clean system for three different temperatures. T\ increases exponentially beyond full polarization 
at the lowest temperature of 50mK. Also note the drop in T\ below full polarization. 



16 




B (Tesla) 



Figure 6: The Polarization P(B) for A = 2.00 and V = WOmK. 

V. QUALITATIVE EFFECT OF LANDAU-LEVEL MIXING 

At fields of the order of a few Tesla the interaction scale e 2 / nl (approximately 
50 VB Kelvin where B is in Tesla) is comparable to the cyclotron scale ftw c (approxi- 
mately 20B Kelvin where B is in Tesla). The ratio scales as l/y/B. The lower the 
field, the more important Landau-level mixing is|33:]. In this section we will show that the 
qualitative effect of Landau-level (LL) mixing is to improve the agreement between theory 
and experiment. 

In previous work we have developed a formalism for taking Landau-level mixing into 
account in our Extended Hamiltonian theory [3^], where the electronic cyclotron coordinate 
rj e is retained. Briefly, by a unitary transformation we eliminate the coupling to the higher 
Landau levels. To leading order, this results in additional four-fermi and six-fermi[3j| terms 

2 

proportional to the perturbation parameter ( = ^J— -• 

While this can be folded into the HF calculation we have described in previous sections, it 
is computationally cumbersome. We will proceed to take its effect into account qualitatively. 
LL-mixing produces level repulsion between the LLL states and higher LL states. The states 
at the top of the LLL are repelled more, with the net effect being to reduce the interaction- 
induced bandwidth of the LLL. This means that the density of states of CF's is increased by 
this effect, implying that T\ should decrease due to LL-mixing. This reduction is larger at 
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low fields and smaller at high fields. Applying this insight to Fig. [3] we see that LL-mixing 
will bring the theoretical prediction into better agreement with the data. 

The magnitude of this effect can be estimated by referring to our previous calculations 
in the gapped fractions[34]. For 3 Tesla, the LL-mi xing parameter is roughly £ = 1.4, while 
for 6 Tesla is decreases to 1. From Fig. 1 of Ref. 34] it is seen that the reduction of the 
FQH gap is between 5-7% in this range. We can thus expect the density of states of the 
CF's to be enhanced a few percent at 3 Tesla relative to the density of states at 6 Tesla. 
This will tend to flatten the concavity of the theoretical curves at low B and improve the 
agreement with experiment. 



VI. CONCLUSIONS, CAVEATS, AND OPEN QUESTIONS 

In this paper we have focussed recent experiments of Tracy and Eisenstein[24] on a /' 
1/2 system for temperatures as low as 45mK and a range of magnetic fields. They measured 
the nuclear relaxation rate 1/T% by magnetically disturbing the system and measuring the 
relaxation of the resistance to its equilibrium value. The predictions of the naive CF picture 
disagree strongly with the measurements. We showed here that the extended Hamiltonian 
theory developed by us gives a satisfactory account of the data if in addition to the thickness 
parameter A of the electron-electron interaction H] we incorporate some constant, spin- 
independent, disorder-induced width of order T = lOOmK for all momentum states. 

Our Hartee-Fock calculation ignores vertex corrections due to gauge fields 

HQ 

since 

they affect properties only in the extreme low-energy/low-temperature limit. We also ig- 
nore vertex corrections due to static disorder because the property we calculate, 1/T 1? is an 
integral over all momenta, and is therefore expected to be insensitive to diffusion-like con- 
tributions at small q. Our optimal value for the disorder width T is lOOmK — 120mK, not 
too different from a crude estimate based on an application of the Drude formula (Eq. (jjj)) 
to the measured longitudinal resistivity Hi- This gives us confidence that our calculation is 
certainly an important part of the explanation. 

The value we have used for the thickness parameter in our calculations, A = 2.00 at 
B = 4 Tesla, requires some discussion. If one took the thickness parameter arising from 



the self-consistent wavefunction in the transverse direction 
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251 ] one would obtain A 



1. The primary effect of A is to soften the Coulomb interaction at short distances. In 
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incompressible states where there is no linear screening, the finite thickness is indeed the 
dominant contribution to softening the Coulomb interaction at short distances. However, in 
a compressible system such as v — 5, there will be further screening by the CF's themselves, 
which will increase the effective value of A beyond that deduced from the finite thickness of 
the 2DEG. We believe this is the reason our value of A is so large. 

Let us now turn to some caveats. The fundamental assumption underlying our approach 
is that disorder averaging gives a true physical picture of the system. While we expect this 
to be true when the system is far from full polarization and the states of both spins are 
extended, we expect that more complicated physics is relevant near full polarization. In this 
regime, we can expect the minority spins to become localized, perhaps in droplets at the 
minima of the disorder potential. Because they are surrounded by a sea of majority spins, in 
addition to the usual disorder potential, the minority spins feel an exchange potential which 
tends to further localize them. Such effects occur in zero-field Fermi liquids as well, but here 
the scale of the exchange interactions is comparable to the Fermi energy, and exchange effects 
are magnified. Thus it is doubtful whether disorder-averaging will give a physically relevant 
answer. We can estimate the behavior of T\ for large B in this regime as follows: Visualize 
droplets of minority spin CF's filling the minima of the disorder potential. As the Zeeman 
energy increases, the droplets will shrink such that their area depends linearly on the Zeeman 
energy. Since the total relaxation rate should be roughly proportional to the total area of 
such droplets (relaxation requires that both spin species of CF's be present at a given nuclear 
site in the localized regime) 1/Ti should decrease linearly with the Zeeman energy. Actually, 
some droplets will disappear when the Zeeman energy increases, so the decrease of 1/Ti will 
be faster than linear in Ez, with the actual power depending on the distribution of minima 
of the disorder potential. Note that the disorder potential which should be used here is not 
the bare one due to inhomogenieties of the dopants, but rather the bare potential screened 
by the compressible CF state, and including exchange contributions of the type mentioned 
earlier. A quantitative calculation of these effects is beyond the scope of this paper. In light 
of these arguments, it is rather fortuitous that our disorder-averaged calculation seems to 
track the data (Fig. [2]) even beyond nominal full polarization. 

Finally, we have made no attempt to construct a theory of the longitudinal 
conductance!^ in the LLL. As has been known for some time 3s|. the dominant terms 
in the current operator are Landau-level mixing terms. However, it seems likely that the 
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low-frequency longitudinal conductance is controlled entirely by LLL physics. An intriguing 
part of the data of TE 2j| is the peak in j—jj^f near nominal full polarization. Our picture 
of the formation of localized droplets of minority spin near nominal full polarization may 
have some relevance for this peak as well. If such droplets are formed, they represent an 
additional repulsive exchange potential for the majority spins, which are mainly responsible 
for transport. Near nominal full polarization, the system is extremely sensitive to the Zee- 
man energy, and we expect the scattering due to the exchange potential of these droplets to 
be maximized. We hope to revisit this issue in future work. 
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